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Abstract 



We present fixed-t subtracted dispersion relations for Compton scattering 
off the nucleon at energies < 500 MeV, as a formalism to extract the 
nucleon polarizabilities with a minimum of model dependence. The subtracted 
dispersion integrals are mainly saturated by vrA^ intermediate states in the 
s-channel — > vrA^ jN and vrvr intermediate states in the t-channel 
77 TTTT NN. For the subprocess 77 — > tttt, we construct a unitarized 
amplitude and find a good description of the available data. We show results 
for Compton scattering using the subtracted dispersion relations and display 
the sensitivity on the scalar polarizability difference a — P and the backward 
spin polarizability 7,^, which enter directly as fit parameters in the present 
formalism. 

PACS : 13.60.Fz, 11.55.Fv, 14.20.Dh, 13.40.-f 

I. INTRODUCTION 

Compton scattering off the nucleon is determined by 6 independent helicity amplitudes 
Ai [i = 1,...,6), which are functions of two variables, e.g. the Lorentz invariant variables u 
(related to the lab energy of the incident photon) and t (related to the momentum transfer 
to the target). In the limit — > 0, the general structure of these amplitudes is governed by 
low energy theorems (LET) based on Lorentz and gauge invariance. These theorems require 
that (I) the leading term in the expansion in u is determined by the global properties of 
the nucleon, i.e. its charge, mass and anomalous magnetic moment, and (II) the internal 
structure shows up only at relative order z/^ and can be parametrized in terms of polariz- 
abilities. In this way there appear 6 polarizabilities for the nucleon, the familiar electric and 
magnetic (scalar) polarizabilities a and /? respectively, and 4 spin (vector) polarizabilities 71 
to 74. These polarizabilities describe the response of the system to an external quasistatic 
electromagnetic field, and as such they are fundamental structure constants of the composite 
system. In particular, these polarizabilities allow one to make contact with classical physics 
phenomena, the dielectric constant and the magnetic permeability of macroscopic media as 
well as the Faraday effect in the case of the spin polarizabilities. 

As a consequence of LET, the differential cross section for u —>■ is given by the (model 
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independent) Thomson term. In a low-energy expansion, the electric and magnetic polariz- 
abilities then appear as interference between the Thomson term and the subleading terms, 
i.e. as contribution of 0(z/^) in the differential cross section, and a and /? can in principle 
be separated by studying the angular distributions. However, it has never been possible to 
isolate this term and thus to determine the polarizabilities in a model independent way. The 
obvious reason is that, for sufficiently small energies, say z/ < 40 MeV, the structure effects 
are extremely small and hence the statistical errors for the polarizabilities large. At larger 
energies, however, the higher terms of the expansion, 0(z/^), become increasingly important. 
Therefore, a reliable theoretical estimate of these higher terms is of utmost importance. 
Moreover, at that order also the spin-dependent polarizabilities come into the game, which 
has the further consequence that a full determination of the 6 polarizabilities will require an 
experimental program with polarized photons and polarized nucleons. 

With the advent of high duty-factor electron accelerators and laser backscattering tech- 
niques, new precision data have been obtained in the 90's and more experiments are expected 
in the near future. In 1991 the Illinois group [|l| measured differential cross sections with 
tagged photons at low energy. As was to be expected, the small counting rate and the 
low sensitivity to structure effects allowed for a reduced statistical precision only. Shortly 
after, Zieger et al. determined the cross section for photon scattering at 6' = 180° by 
detecting the recoil proton in a magnetic spectrometer. In a series of experiments, the 
Illinois-Saskatoon group then studied angular and energy distributions over a wider range. 
Hallin et al. investigated the region from pion production threshold to the A(1232) reso- 
nance with a high duty-factor bremsstrahlung beam. Though the statistical and systemat- 
ical errors were small, the range of energy was clearly outside of the low-energy expansion. 
The presently most accurate values for the proton polarizabilities were derived from the 
work of MacGibbon et al. whose experiments were performed with tagged photons at 
70 MeV< V < 100 MeV and untagged ones at the higher energies, and analyzed in collabo- 
ration with L'vov 1^ by means of dispersion relations (in the following denoted by DR) at 
constant t. The final results were 

a = (12.1 ± 0.8 ± 0.5) X 10"^ fm^ , 

p = (2.1 ^ 0.8 T 0.5) X 10"^ fm^ . (1) 



The physics of the A(1232) and higher resonances has been the objective of further re- 
cent investigations with tagged photons at Mainz and with laser-backscattered photons 
at Brookhaven The measured differential cross sections and polarization asymmetries 
helped to discard a long-standing problem, a unitarity violation reported in earlier exper- 
iments, and provided useful cross-checks for the magnetic dipole and electric quadrupole 
excitation of the A resonance. Such data were also used to give a first prediction for 
the so-called backward spin polarizability of the proton, i.e. the particular combination 
Ttt = 7i + 72 + 274 entering the Compton spin-flip amplitude at 6* = 180° [||, 



7^ 



_j_2 g 

27.1 ± 2.2(stat + syst) ^'^ (model) 



X 10^^ fm' 



(2) 



In 1991 Bernard et al. |^ evaluated the one-loop contributions to the polarizabilities in the 
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framework of relativistic chiral perturbation theory (ChPT), with the result a= 10-/3 = 12.1, 
(here and in the following, the scalar polarizabilities are given in units of 10"'* fm^). In or- 
der to have a systematic chiral power counting, the calculation was then repeated in heavy 
baryon ChPT |10| to O(p^), the expansion parameter being an external momentum or the 



quark mass. The result explained the small value of the magnetic polarizability, which had 
been difficult to obtain in quark model calculations. A further calculation to 0{j)^^ resulted 
in the values a = 10.5 ±2.0 and (3 = 3.5 ±3.6, the errors being due to 4 counter terms enter- 



ing to that order, which were estimated by resonance saturation |Tl|. One of these counter 
terms describes the paramagnetic contribution of the A(1232), which is partly cancelled by 
large diamagnetic contributions of pion-nucleon loops. 



In view of the importance of the A resonance, Hemmert, Holstein and Kambor ||T^ pro- 
posed to include the A as a dynamical degree of freedom. This added a further expansion 
parameter, the difference of the A and nucleon masses ("e expansion"). A calculation to 
0{e^) yielded the results (see also Ref. 0) 

a = 12.2 ±0 ±4.2 = 16.4 , 

/5 = 1.2 ±7.2 ±0.7 = 9.1 , (3) 

where the 3 terms on the rhs are the contributions of pion-nucleon loops (identical to the 
predictions of the 0{p^) calculation), A pole terms, and pion-A loops. These predictions 
are clearly at variance with the data, in particular a + (3 = 25.5 is nearly twice the rather 
precise value determined from DR (see below). In an optimistic view the problem is to be 
cured by an O(e^) calculation, which is likely to produce a large diamagnetism as observed 
by Ref. [0. On the other hand, a pessimist might doubt that the expansion converges 
sufficiently well before the higher orders have introduced a host of unknown counter terms. 



The spin polarizabilities have been calculated in both relativistic one-loop ChPT 1 10,15 



and heavy baryon ChPT |]T4|. In the latter approach the predictions are 



70 = 4.6 -2.4- 0.2 ±0 = ±2.0 , 

7^ = 4.6 ±2.4 -0.2 -43.5 = -36.7 , (4) 

where the spin polarizabilities are given here and in all of the following in units of 10~^ fm^. 
The 4 separate contributions on the rhs of Eq. (^ refer to Nvr-loops, A-poles, A7r-loops, and 
the triangle anomaly, in that order. It is obvious that the anomaly or 7r°-pole contribution 
is by far the most important one, and that it would require surprisingly large higher order 
contributions O(e^) to increase to the value of Ref. M. Similar conclusions were reached 



in the framework of DR. Using the framework of DR at t = const of Ref. fg], Ref. ||T6 



obtained a value of 7,^ = —34.3, while L'vov and Nathan [l?! worked in the framework of 
backward DR and predicted 7^ = —39.5 ±2.4. In the latter approach the dispersion integral 
is drawn along a line t(z/) corresponding to backward Compton scattering, i.e. on the lower 
boundary of the physical s-channel region, which is then complemented by a path into the 
physical t-channel region. 

As we have stated before, the most quantitative analysis of the experimental data has been 
provided by dispersion relations. In this way it has been possible to reconstruct the forward 
non spin-flip amplitude directly from the total photoabsorption cross section by Baldin's 
sum rule |lH], which yields a precise value for the sum of the scalar polarizabilities 
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a + (3 = 14.2 ± 0.5 (Ref. [19]) 

= 13.69 ± 0.14 (Ref. |2|) . 



(5) 



Similarly, the forward spin-flip amplitude can be evaluated by an integral over the difference 
of the absorption cross sections in states with helicity 3/2 and 1/2, 



While these predictions rely on pion photoproduction multipoles, the helicity cross sections 
have now been directly determined by scattering photons with circular polarizations on po- 
larized protons ||22|| . 

In the case of forward Compton scattering the momentum transfer and, hence, the Man- 
delstam variable t vanishes. In that sense the above sum rules of Eqs. (^, ^ are derived 
from DR at t = 0. At finite angles, however, the analysis requires DR at t = const. < 0, in a 
range of values between (forward scattering) and the largest negative value oit {t = tmax) , 
determined by the largest scattering angle at the highest photon energy. As mentioned 
above, the most quantitative and detailed such analysis has been performed by L'vov and 
collaborators |^,^ in the framework of unsubtracted DR at t = const. Unfortunately, not 
all of the dispersion integrals converge, as can be inferred from Regge theory. The reason 
for the divergence of the integrals is essentially given by fixed poles in the t channel, notably 
the exchange of the neutral pion and of a somewhat fictitious a meson with a mass of about 
600 MeV and a large width, which models the two-pion continuum with the quantum num- 
bers / = J = 0. In a more formal view, the dispersion integral is performed along the real 
axis in the range —Vmax < ^ < +Vmaxi with v^ax ~ 1-5 GeV, and then closed by a semi- 
circle with radius Vmax in the upper half of the complex i/-plane. The contribution of the 
semi-circle is then identified with the asymptotic contribution described by t-channel poles. 
This introduces unknown vertex functions and the mass of the "a meson" , which have to be 
determined from the experiment. Moreover, the analysis depends appreciably on the choice 
of I'maxi and there are substantial contributions of intermediate states beyond the relatively 
well-known pion-nucleon continuum. These higher states include multipion, rj- and p-meson 
production, Avr-loops and nonresonant s-wave background. The physics behind these effects 
is certainly worthwhile studying, and there can be no doubt that within the next years we 
shall learn more about them by detailed coincidence studies of multipion and heavier meson 
production, but also directly from a careful analysis of Compton scattering at the higher 
energies 0. However, the quest for the polarizabilities as fundamental structure constants 
should not be burdened by too many open questions and phenomenological models. 

In view of the problems of unsubtracted DR, we propose to analyse Compton scattering 
in the framework of subtracted DR at constant t, with the eventual goal to determine the 
6 polarizabilities with the least possible model dependence. We choose to subtract the 6 
Compton amplitudes Aj(z/, t) at the unphysical value u = 0, i.e. write subtracted DR for 
y4j(z/, t) — Ai{0, t) at constant t. As we shall show in the following, these subtracted DR con- 
verge nicely and are quite well saturated already at ~ 400 — 550 MeV, i.e. essentially by 
one-pion production. Clearly the price to pay are 6 new functions Ai{0, t), which have to be 
determined by another set of dispersion relations, at z/=const=0 and by use of information 
obtained from the t-channel reaction 77 anything. 



7o = 7i - 72 - 274 



1.34 (Ref. [21]) 
0.6 (Ref. [0) . 



(6) 
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In order to reduce the dependence on the higher intermediate states in the t-channel, we 
subtract again, i.e. write DR for Aj(0, t) — Aj(0, 0), the subtraction constants Aj(0, 0) being 
hnear combinations of the 6 polarizabihties. Since 4 of these subtraction constants can be 
calculated from unsubtracted DR at )f:=const, only 2 parameters have to be fixed by a fit to 
low energy Compton scattering, the combinations a — P and 7,1 = 7i + 72 + 274 describing 
the backward non spin-flip and spin-flip amplitudes, respectively. 

In a somewhat similar approach, Holstein and Nathan [^] combined s- and t-channel in- 
formation to predict the backward scalar polarizability a — (3. Using unsubtracted backward 
DR they obtained, from the integration along the lower boundary of the s— channel region, 
the result [a — (3Y = —6 ± 3, and from the t-channel region a contribution of {a — jSy ~ 9. 
The sum of these two contributions, a — /3 ~ 3 ±3, is at variance with the presently accepted 
experimental (global average) value, a — P = 10.0 ± 1.5 ± 0.9 of Ref. [Q. The difficulty to 
predict this observable is due to the bad convergence of the integrals in both the s— and the 
t-channel regions. As may be seen from Fig. 3 of Ref. [^, the t-channel integral obtains 
quite sizeable contributions from 10 < t < 40 m^, in which region the integrand changes 
its sign. Independent of the numerical analysis, the authors find an extremely interesting 
relationship connecting the polarizabihties of pions and nucleons. This connection comes via 
the t-channel integral for the nucleon and the low-energy expansion of the s-wave 77 — *• tttt 
amplitude, and results in 6a = —5(3 = 0.5a'^ ~ 1.4, which is the contribution of the pion 
polarizability a'^ to the nucleon's electric and magnetic polarizabihties. In concluding this 
discussion we point out the difference of our present approach and the calculation of Ref. p4 



We do not intend to predict a — p. Instead we want to develop a dispersion description of 
Compton scattering allowing for a derivation of a — /? with a minimum of model dependence. 
For this purpose we use a scheme of subtracted DR whose subtraction constants are linear 
combinations of polarizabihties, which have to be determined by a fit to the Compton data. 

In section |I| we shall give a general introduction to subtracted DR. This technique is then 
applied to the cases of DR at t=const (s-channel dispersion integral) and DR at u = 
(t-channel dispersion integral) in sections |T| and respectively. Our results are compared 
to the existing low-energy Compton data in section ^ and our conclusions are drawn in 
section [V^. 



II. FIXED-T SUBTRACTED DISPERSION RELATIONS 

Assuming invariance under parity, charge conjugation and time reversal symmetry, the 
general amplitude for Compton scattering can be expressed in terms of six independent 
structure functions Ai{i/,t), i = I, ...,6. These structure functions depend on two Lorentz 
invariant variables, e.g. ly and t as defined in the following. Denoting the momenta of the 
initial state photon and proton by q and p respectively, and with corresponding final state 
momenta q' and p', the familiar Mandelstam variables are 

s = {q + p)\ t = {q-q')\ u = {q-p'f. (7) 

These variables fulfill the constraint 

s + t + u = 2M^ . (8) 
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The variable v is defined by, 



V = = E^^ , (9) 

where is the photon energy in the lab frame and M the nucleon mass. The Mandelstam 
plane is shown in Fig. 1, and the boundaries of the physical and spectral regions are dis- 
cussed in Appendix A. 

The invariant amplitudes Ai are free of kinematical singularities and constraints, and 
because of the crossing symmetry they satisfy the relation y4i(z/, t) = Ai{—u,t). Assum- 
ing further analyticity and an appropriate high-energy behavior, the amplitudes fulfill 
unsubtracted dispersion relations at fixed t, 

ReAiu,t) = Afiu,t) + - V dv' ' , 10 

where Af are the Born (nucleon pole) contributions, Im^Aj the discontinuities across the 
s-channel cuts of the Compton process and vthr = ''^tt + {fnl + t/2)/(2M). However, such 
unsubtracted dispersion relations require that at high energies (z/ oo) the amplitudes 
lmsAi{i',t) drop fast enough so that the integral of Eq. (|TUp is convergent and the contribu- 
tion from the semi-circle at infinity can be neglected. For real Compton scattering, Regge 
theory predicts the following high-energy behavior for z/ — > oo and fixed t 0: 



y4i 2 ~ z/' 



a(i) 
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^,5,6 ~ , A, ~ Z/-W-3 , (11) 

where a(t) < 1 is the Regge trajectory. In particular we note that the Regge trajectory 
with the highest intercept, i.e. a(0) ~ 1.08, corresponds to soft pomeron exchange. Due 
to this high energy behavior, the unsubtracted dispersion integral of Eq. (|T0|) diverges for 
the amplitudes Ai and A2. In order to obtain useful results for these two amplitudes, L'vov 
et al. proposed to close the contour of the integral in Eq. (^) by a semi-circle of finite 
radius Umax (instead of the usually assumed infinite radius!) in the complex plane, i.e. the 
real parts of Ai and A2 are calculated from the decomposition 

ReA{u,t) = (z/,t) + A^iu.t) + Ariiy,t) , (12) 

with A™* the s-channel integral from pion threshold Uthr to a finite upper limit Umax, 

AfM^ivr,.'"-^^^^, (13) 



and an 'asymptotic contribution' representing the contribution along the finite semi- 
circle of radius Vmax in the complex plane. In the actual calculations, the s-channel integral 
is typically evaluated up to a maximum photon energy E^ = Vmaxif) — t/{4:M) ^ 1.5 GeV, 
for which the imaginary parts of the amplitudes can be expressed through unitarity by 
the meson photoproduction amplitudes (mainly Ivr and 2n photoproduction) taken from 
experiment. All contributions from higher energies are then absorbed in the asymptotic 
term, which is replaced by a finite number of energy independent poles in the t-channel. In 
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particular the asymptotic part of Ai is parametrized in Ref. |^ by the exchange of a scalar 
particle in the t-channel, i.e. an effective "cr meson" , 



ATM^AUt) = EzTl9^^ (14) 

t — mi 



where rria- is the "a meson" mass, and g^NN and F^.^^ are the couplings of the "a meson" to 
the nucleons and photons respectively. The asymptotic part of A2 is parametrized by the 
vr^ t-channel pole. 

This procedure is relatively save for A2 because of the dominance of the vr" pole or tri- 
angle anomaly, which is well established both experimentally and on general grounds as 
Wess-Zumino-Witten term. However, it introduces a considerable model-dependence in the 
case of Ai. Though "cr mesons" have been repeatedly reported in the past, their prop- 
erties were never clearly established. Therefore, this particle should be interpreted as a 
parametrization of the I = J = part of the two-pion spectrum, which shows up differently 
in different experiments and hence has been reported with varying masses and widths. 

It is therefore the aim of our present contribution to avoid the convergence problem of 
unsubtracted DR and the phenomenology necessary to determine the asymptotic contribu- 
tion. The alternative we shall pursue in the following is to consider DR at fixed t that are 
once subtracted at u = 0, 

ReA{u,t) = Af{u,t) + [A{0,t) - Af{0,t)] + - V du' f,,'^ ';! • (15) 

These subtracted DR should converge for all six invariant amplitudes due to the two ad- 
ditional powers of u' in the denominator, and they are essentially saturated by the vrA^ 
intermediate states as will be shown in section |T|. In other words, the lesser known con- 
tributions of two and more pions as well as higher continua are small and may be treated 
reliably by simple models. 

The price to pay for this alternative is the appearance of the subtraction functions 
Ailiy = 0,t), which have to be determined at some small (negative) value of t. We do 
this by setting up a once-subtracted DR, this time in the variable t, 

A{0,t) - Af{0,t)=[A{0,0) - Af(0,0)] + [A'rP"'%0,t) - A*-^°'^(0,0) 

t_ Im,A.(0,tO t_ r Im,A.(0,tO 

7(2^^)2 t'{t'-t) TT t'{t'-t) ' ^ ' 

where A-~^°''^(0, t) represents the contribution of poles in the t-channel, in particular of the 
TT^ pole in the case of A2, which is given by 



Af(0,t) = (17) 



The coupling determined through the tt^ 77 decay as 



= ^rnloF^,. (18) 
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Using r (vr^ — > 77) = 7.74 eV one obtains -F^o^-y = -0.0252 GeV ^, where the sign is 
in accordance with the 7r°77 couphng in the chiral hmit, given by the Wess-Zumino-Witten 



effective chiral Lagrangian. The irNN couphng is taken from Ref. |2^ : g^^j^/Air = 13.72. 
This yields then for the product of the couplings in Eq. (|l^) : F^o^^ g^^^j^ pa -0.331 GeV~^. 

The imaginary part in the integral from 4m^ — > 00 in Eq. (^) is saturated by the possible 
intermediate states for the t-channel process (see Fig. |]), which lead to cuts along the positive 
t-axis. For values of t below the KK threshold, the t-channel discontinuity is dominated 
by TTTT intermediate states. The second integral in Eq. (|16D extends from —00 to a, where 
a = —4 (m^ + 2Mm^) ^ —1.1 GeV^ is the boundary of the su spectral region for i/ = 
(see Appendix ^ for a detailed discussion). As we are interested in evaluating Eq. (|T6|) for 
small (negative) values of t (|t| <^ \a\), the integral from —00 to a will be highly suppressed 
by the denominator of the subtracted DR, and can therefore be neglected. Consequently, 
we shall saturate the subtracted dispersion integrals of Eq. ([I6| ) by the contribution of vrvr 
intermediate states, which turns out to be a good approximation for small t. We will show 
the convergence of the t-channel dispersion integral in section |^ and thus verify the quality 
of the approximation. 

The t-dependence of the subtraction functions Ai{0,t) is now determined, and only the 
subtraction constants y4j(0,0) remain to be fixed. We note that the quantities 

a, = Ai{0,0) - Af(0,0) (19) 

are directly related to the polarizabilities, which can then be obtained from a fit to the 
Compton scattering data. For the spin-independent (scalar) polarizabilities a and /3, one 
finds the two combinations 

a + (3 = [as + ae) , (20) 

= -7^«i, (21) 

which can be determined from forward and backward scattering respectively. Furthermore, 
the forward combination a + P is related to the total absorption spectrum through Baldin's 
sum rule |T^, 



00 



, a{^N ^ X) 



{a + PU = ^ / du'^-^. (22) 

" ^thr 



The 4 spin dependent polarizabilities 71 to 74 of Ragusa |^ are defined by 



7o = 7i - 72 - 274 = 04 , (23) 

713 = 71 + 273 = - (as + ae) , (24) 

714 = li- 274 = (2 a4 + as - ae) , (25) 
In = 71 + 72 + 274 = - (02 + as) , (26) 

where 70 and 77r are the spin (vector) polarizabilities in the forward and backward directions 
respectively. Since the 7r° pole contributes to A2 only, the combinations 70, 713 and 714 of 
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Eqs. (p3D-(p5D are independent of this pole term [jl6[, and only the backward spin polariz- 
ability 7^ is affected by the anomaly. 

Although all 6 subtraction constants ai to of Eq. (|1^) could be used as fit parameters, 
we shall restrict the fit to the parameters ai and 02, or equivalently to a — 13 and 7^. The 
subtraction constants 04, as and qq will be calculated through an unsubtracted sum rule, as 
derived from Eq. (p!0|) . 



04,5,6 



2 



Im,A4,5,6(^^t=0) 



(27) 



The remaining subtraction constant 03, which is related to a + f3 through Eq. 
fixed through Baldin's sum rule, Eq. 
13.69. 



will be 



using the value obtained in Ref. : a + /3 



III. S-CHANNEL DISPERSION INTEGRAL 

In this section we describe the calculation of the s-channel contributions, which enter in 
the once-subtracted dispersion integral of Eq. (|T5|) and in the calculation of subtraction 
constants 04, 05 and ag through Eq. (^). The imaginary part of the Compton amplitude 
due to the s-channel cuts is determined from the scattering amplitudes of photoproduction 
on the nucleon by the unitarity relation 

2Im,Ty, = 5^(27r)V(Px - Pi)TlfTxi , (28) 

X 

where the sum runs over all possible states that can be formed in the photon-nucleon reac- 
tion. Due to the energy denominator l/i/'(z/'^ — z/^) in the subtracted dispersion integrals, 
the most important contribution is from the vrA^ intermediate states, while mechanisms 
involving more pions or heavier mesons in the intermediate states are largely suppressed. 
In our calculation, we evaluate the vrA^ contribution using the multipole amplitudes from 
the analysis of Hanstein, Drechsel and Tiator p8[ at energies E-^ < 500 MeV and at the 
higher energies we take as input the SAID multipoles (SP98K solution) The expan- 

sion of lm.sAi into this set of multipoles is truncated at a maximum angular momentum 
jmax = / ± 1/2 = 7/2, with the exception of the lower energy range {E^ < 400 MeV) where 
we use jmax = 3/2. The higher partial waves with j > jmax + 1 are evaluated analytically 
in the one-pion exchange (OPE) approximation. The relevant formulas to implement the 
calculation are reported in Appendix B and C of Ref. . 

The multipion intermediate states are approximated by the inelastic decay channels of the 
T^N resonances. In the spirit of Ref. |^ and the more recent work of Ref. we assume 



that this inelastic contribution follows the helicity structure of the one-pion photoproduction 
amplitudes. In this approximation, we first calculate the resonant part of the pion photo- 
production multipoles using the Breit-Wigner parametrization of Ref. [^], which is then 



scaled by a suitable factor to include the inelastic decays of the resonances. The resulting 
contribution to Im^ Ai is 

[Im, = R [ini, Aif*-''''^ , (29) 



9 



with the ratio R given by 



In Eq. (|30D , B^^ is the single-pion branching ratio of the resonance N* and r^(iy) the 
energy-dependent pionic width while the inelastic width finei(H^) of the decays 

N* —>■ {tttcN, rjN, tctcttN,..) is parametrized as in Ref. in order to provide the correct 
threshold behavior for the resonant two pion contribution. 

The ttN channel consistently reproduces the measured photoabsorption cross section in the 
energy range < 500 MeV, while at the higher energies nonresonant mechanisms should 
be included in addition to the resonant mechanism to fully describe the multipion channels. 
In Ref. P], the non- resonant contribution to the two-pion photoproduction channel was ap- 
proximately taken into account by calculating the OPE diagram of the — > vrA reaction. 
The difference between the data and the model for two-pion photoproduction consisting of 
resonant mechanisms plus the OPE diagram for the nonresonant mechanism, was then fitted 
in Ref. and attributed to a phenomenological, non-resonant —>■ vrA s-wave correction 
term. 

A more detailed description of the vrvrA^ channel is clearly worthwile to be undertaken, 
especially in view of the new two-pion photoproduction data (both unpolarized and polar- 
ized) that will be available from MAMI and JLab (CLAS) in the near future. However, 
for the extraction of the polarizabilities, the strategy followed in this paper is to minimize 
sensitivity and hence model uncertainty to these higher channels. 

We show in Fig. |] that in the subtracted DR, the sensitivity to the multipion channels 
is indeed very small. For the unsubtracted DR, on the other hand, the influence of the 
multipion channels amounts to about 30 % of the amplitude A2. We furthermore note from 
Fig. ^that the subtracted DR are essentially saturated at ~ 0.4 GeV and that they only 
receive a negligible contribution from multipion channels. The importance of the multipion 
channels is even weaker in the case of the amplitudes ^3 to Aq. 

In Table I and II, we show our predictions for the dispersion integral of the spin polar- 
izabilities of the proton and neutron, respectively. We list the separate contribution of the 
ttN channel, HDTflvr), and the total result, HDT, which includes the inelastic resonance 



channels. The last column shows the values of the dispersion calculation of Ref. , which 
is based on the one-pion multipoles of the SAID-SP97K solution and the model for double- 
pion production mentioned above. The small differences between the one-pion multipoles of 
SAID-SP97K and SAID-SP98K at the higher energies are practically negligible for the spin 
polarizabilities, while the results are very sensitive to the differences between the HDT and 
SAID analyses. As discussed in Ref. this fact is mainly due to a different behaviour of 
the i?o+ partial wave near threshold, giving rise to substantial effects in the case of the for- 
ward spin polarizability. While the one-pion contribution from SAID-SP98K is 79 = —1.26 
and 7o = —0.03, we obtain 7q = —0.75 and 79 = —0.06 with the HDT multipoles for 
< 500 MeV. 
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IV. T-CHANNEL DISPERSION INTEGRAL 



We next evaluate the t-channel dispersion integral in Eq. (|T6D from 4m^ to oo. The 
kinematics of the t-channel reaction 77 NN is shown in Fig. ^ The subtracted dispersion 
integral is essentially saturated by the imaginary part of the t-channel amplitude 77 — ^ NN 
due to vrvr intermediate states. To calculate this contribution, we have to construct the 
amplitudes 77 — *• tttt and vrvr — » NN. 

We start with the isospin and helicity structure of the 77 —>■ tttt amplitude, denoted by F. 
Because of the Bose symmetry of the 77 state, only the even isospin values 1 = and 2 are 
possible. We can express the charged (77 tt+tt^) and neutral (77 tt^tt^) amplitudes in 
terms of those with good isospin by 

= + (charged pions) , 

= -^1^'=° + \/f^'^' (neutral pions) . (31) 

The reaction 77 nn has two independent helicity amplitudes Fj>^^{t, 6 t^t^), where = 
A^^ — A^, being the difference of the final photon helicity {X'^) and the initial photon helicity 
(A^), takes on the values or 2, depending upon whether the photons have the same (A^ = 
0) or opposite (A^ = 2) helicities. The 77 — >• vrvr helicity amplitudes depend upon the cm. 
energy squared t, and the pion cm. scattering angle ^^^r- In terms of the helicity amplitudes 
F\^, the 77 — s> TTTT differential cm. cross section is given by 

da \ P 



d cos 9^ 



) = {l^A,=o(t,^..)P + |FA,=2(t,^..)P} , (32) 

/cm. 647rt 



with /3 = a/1 — 4m^ / 1 the pion velocity. In Appendix 0, we give the partial wave expansion 
of the 77 — >• TTTT helicity amplitudes Fl_^(t, 6*7^77) for a state of isospin J, and thus define the 
partial wave amplitudes Fj^ (t) (see Eqs. ( [BToD and (pTsD), where J can only take on even 
values. 

To construct the helicity amplitudes Fa^ for the process 77 tttt, we first evaluate 
the Born graphs as shown in Fig. |. These graphs only contribute to the charged channel 
77 TT+TT^. The Born contributions to the helicity amplitudes F^;^ ^ are denoted as B\ 
and given by 

1 /52 
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B 



(33) 

The partial wave expansion of the Born terms BjA^it) is discussed in Appendix ^ 
(Eq. ( |B16| )). As the Born amplitudes are only non-zero for the charged pion channel, the 
two isospin amplitudes of Eq.(0) are related by 
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We now construct the unitarized amplitudes Fjj^^{t), starting from the Born amphtudes 
(t) and following the method outlined in Refs. [0,0 • We first note that the low energy 
theorem requires for each partial wave that 



B 



1 , as t 



0. 



(35) 



Next, the invariant amplitude for the process 77 vrvr is assumed to have Mandelstam 
analyticity. Each partial wave then has a right-hand cut from t = 4m^ to +00 and a 
left-hand cut from t = —00 to 0. Though the Born amplitude is real for all values of t, 
its partial waves are complex below t = 0. The partial waves of the full amplitude have no 
other sources of complexity in this region, and so we can write a DR for the difference of 
the full and the Born amplitudes. 



JA. 



t{t - Ami 



dt' 

Ami t'{t' 



ImFj, (f) 



Aml)^{t' -t) 



(36) 



with an additional factor of {t{t — Am'^)i)~^ providing the right asymptotics for the con- 
vergence of the integral. The next step is to evaluate the imaginary part of the amplitude 
in Eq. (|3^) . To do this, we exploit the unitarity condition 



ImFj^^(77 ^ tttt) = 5]]p„FX(77 



n) l]{n 



TTTT] 



(37) 



where pn are the appropriate kinematical and isospin factors for the intermediate channels 
n, and X(n — tttt) is a hadronic amplitude. Below the next inelastic threshold, it follows 
from unitarity that the phase iiJ^'^'^) r-f +1^^ -r^inoc^ 

of the corresponding tttt tttt partial wave. 



yj of each partial wave Fj^ is equal to the phase S^.^ 



ImFj, (77 



TTTT J = p 
^/ (77 tttt) 



tttt) Xj(nn 



nn 



(38) 



This fact can be incorporated into the so-called Omnes function, which is constructed to 
have the phase of the vrvr scattering amplitude above vrvr threshold, and to be real otherwise, 



exp 



dt' —— 

Ami ^ 



IE] 



(39) 



The function Fjy^^(f2j) ^{t) is by construction real above vrvr threshold, but complex below 
threshold due to the complexity of the Born partial waves Bj^^ . Hence we can write a 



dispersion relation for 



^ JA^ 



BjA^ 



{n'jr\t)/t{t-4miy/\ 



t{t - Ami 



2\J/2 



vr 



dt 



Ami 



, B^j^^{t')lm [{n^j)-\t')] 
t'{t' - Amiyi\t' - t) 



■ (40) 
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For t > 4:m1, this integral is understood to be a principal value integral, which we implement 
by subtracting the integrand at t' = t. In this way we obtain a regular integral, which can 
be performed without numerical problems. 



t 



TT 



dt' 
t'{t' - t) 



{t - 4m2)-^/2 



■ (41) 



In our formalism, the s( J = 0)- and d( J 



phaseshifts, we use the solutions that were determined in Ref. |35 



2)-waves are unitarized. For the s- and d-wave tttt 

For the higher partial 



waves, the corresponding rni phaseshifts are rather small and are not known with good 
precision. Therefore, we will approximate all higher partial waves (J > 4) by their Born 
contribution. The full amplitudes for the charged and neutral channels can then be cast 
into the forms 



;i(7r+7r ) 



J=0,2 



J2 V^ITii 



J=0,2 



(J + A^)! 



7/=0 
JA-i 



Xit) - BjA,{t) 
2 



PpicosO^^) , (42) 



it) + 



it) 



P7(cos0, 



(43) 



The two-pion intermediate contribution holds to good precision up to KK threshold 
(~ 1 GeV^), because the four-pion intermediate state couples only weakly and gives only 
small inelasticities in the tttt phaseshifts. 

In Figs. § and 0, we show our results for the total and differential 77 ti^ti~ cross sec- 
tions and a comparison to the existing data. In the threshold region, the charged pion cross 
sections are clearly dominated by the Born graphs of Fig. ^ because of the vicinity of the 
pion pole in the t-channel of the 77 — > tt+tt^ process. However, the results for the unitarized 
calculation show that s-wave rescattering is not negligible but leads to a considerable en- 
hancement at energies just above threshold. Besides the low energy structure, driven by the 
Born terms, the 77 — >■ tttt process has a prominent resonance structure at higher energies 
corresponding to excitation of the isoscalar /2(1270) resonance, with mass nif^ = 1275 MeV 
and width Tf^ = 185.5 MeV p5[. The /2 resonance shows up in the partial wave Fj=2A^=2 
as outlined in Appendix p. Therefore, the most efficient way to unitarize this particular 
partial wave is to make a Breit-Wigner ansatz for the /2 excitation, which is described in 
Appendix P where we also give some details of the formalism for a spin-2 particle. The 
Breit-Wigner ansatz for the /2 contribution to the partial wave Fj=2K~,=2 depends upon the 
couplings /2'7r7r and /277- The coupling /2vr7r is known from the decay /2 tttt and is 
taken from Ref. [^. The coupling /277 is then fitted to the 77 tttt cross section at the 



/2 resonance position, and is consistent with the value quoted in Ref. ||2^. The resulting 
amplitude, consisting of unitarized s-wave, /2 excitation and Born terms for all other partial 
waves (with J > 4) is seen from Figs. ^ and to give a rather good description of the 
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77 vr+7r~ data up to Wt,t, ^ 1.8 GeV. Only in the region Wmr ~ 0.7 - 0.8 GeV, does our 
description slightly overestimate the data. 

Having constructed the 77 — > tttt amplitudes, we next need the tttt — > NN amplitudes 
in order to estimate the contribution of the vrvr states to the t-channel dispersion integral 
for Compton scattering. As we only kept s- and d- waves for 77 tttt, we will only need 
the s- and d-waves (J = 0, 2) for tttt —>■ NN. For every partial wave J, there are two 
independent tttt —>■ NN helicity amplitudes depending on whether the nucleon and 

anti-nucleon have the same {f:l{t)) or opposite {f:[{t)) helicities. We refer the reader to 
Appendix |B| (Eqs. ( |B11 ) and ( |Bl^j| )) for details. In this work, we take the s- and d-waves 



from the work of Hohler and collaborators |3^, in which the lowest tttt NN partial wave 
amplitudes were constructed from a partial wave solution of pion-nucleon scattering, by use 
of the TTTT phaseshifts of Ref . , which we also used to construct the 77 tttt amplitudes. 
In Ref. the tttt NN amplitudes are given for t values up to t ~ 40 ■ ^0.78 GeV^, 
which will serve well for our purpose since the subtracted t-channel dispersion integral will 
have converged much below this value as shown in the following. 

Finally, we can now combine the 77 tttt and tttt NN amplitudes to construct the 
discontinuities of the Compton amplitudes across the t-channel cut. In Appendix 0, we show 
in detail how the Compton invariant amplitudes Ai,...,Aq are expressed by the t-channel 
(77 — > NN) helicity amplitudes. Through unitarity we then express the imaginary parts of 
these t-channel (77 — > NN) helicity amplitudes in terms of the 77 tttt and tttt —>■ NN 
amplitudes. We finally express the discontinuities IratAi of the invariant amplitudes Ai {i 
= 1,...,6) in terms of the corresponding 77 tttt and vrvr NN partial wave amplitudes 
(see Eq. ( |B17| )). As we restrict ourselves to s- and d-wave intermediate states in the actual 
calculations, we give here the expressions ai u = 0, including s- and d-waves only, that are 
needed for the subtracted t-channel dispersion integral of Eq. (^) , 



ImtA2(i/ = 0,t)'" = 0, 
\mtAi{u = 0,t)2^ = , 

lm,A,{u = 0,t)'- = - ( ^Zi^y^' M =o(t) f-{t) 



t V 2 



lm,A,{u = 0,t)'- = - ( ^I±^Y' M ^ F2A^=2(t) fJ{t) . (44) 



t J 2 

The reader should note that the s-wave tttt intermediate state only contributes to the am- 
plitude Ai. It is the t-dependence of this I = J = nir state in the t-channel that is 
approximated in Ref. ^ and parametrized by a "sigma" pole. The d-wave tttt intermediate 
state gives rise to imaginary parts for the amphtudes 741,^43, A^ and Aq. The amplitude A2 
(at z/ = 0) corresponds to the t-channel exchange of an object with the quantum numbers 
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of the pion (e.g. vr" pole in Eq. (|T^). Therefore two-pion intermediate states do not have 
the quantum numbers to contribute to the amphtude A2. The imaginary part of A4 receives 
only contributions from vrvr intermediate states with J > 4 (see Eq. ( [B17| )) and therefore is 
zero in our description, as we keep only s- and d-waves. 

In Fig. ^we show the convergence of the t-channel integral from 4m^ to 00 in the subtracted 
DR of Eq. (|16D . We do so by calculating the dispersion integral as function of the upper 
integration limit tupper and by showing the ratio with the integral for tupper = 0.78 GeV^. 
The latter value corresponds to the highest t value for which the vrvr — > NN amplitudes are 
given in Ref. One clearly sees from Fig. |^ that the unsubtracted t-channel DR shows 
only a slow convergence, whereas the subtracted t-channel DR has already reached its final 
value, within the percent level, at a t value as low as 0.4 GeV^. 



V. RESULTS AND DISCUSSION 

In this section we shall present our results for Compton scattering off the nucleon in the 
dispersion formalism presented above. 

The real and imaginary parts of the six Compton amplitudes are displayed in Fig. |^. 
Note that for the real part, we only show the subtracted s-channel integral of Eq. (P^. 
As can be seen from Fig. |^, these amplitudes show strong oscillations due to interference 
effects between different pion photoproduction multipoles, in particular for threshold pion 
production by -Eo+ and A-excitation by Mi+. 



In Figs. |10| and |TT| we show our predictions in the subtracted DR formalism and compare 
them with the available Compton data on the proton below pion threshold. These data were 
used in Ref. to determine the scalar polarizabilities a and (3 through a global fit, with 
the results given in Eq. (|I]). In the analysis of Ref. |Q, the unsubtracted DR formalism was 
used and the asymptotic contributions (Eq. ([T2|)) to the invariant Compton amplitudes Ai 
and A2 were parametrized. In particular, Ag'^ was described by the 7r° pole, which yields 
the value — —45. The free parameter entering in AI'' was related to a — f3, for which the 
fit obtained the value a — /? ~ 10. Keeping a — (3 fixed at that value, we demonstrate in 



Fig. |T0| that the sensitivity to 7,^ is not at all negligible, especially at the backward angles 
and the higher energies. Although we do not intend to give a best fit at the present stage, the 
subtracted DR formalism allows one to directly use the values a — (3 and 7,^ as fit parameters. 



as is obvious from Eqs. (15, 16). We investigate this further in Fig. 11, where we show our 
results for different a — (3 and for a fixed value of 'Jt^ = —37, which is consistent with the 
heavy baryon ChPT prediction [|l^ and close to the value obtained in Ref. |jl^ in a backward 



DR formalism. For that value of 77r, a better description of the data (in particular at the 
backward angle) seems to be possible by using a smaller value for a — f3 than determined 
in Ref. [Q. For a more reliable extraction of the polarizabilities, more accurate data over 
the whole angular and energy range are necessary. Recently, Compton data were taken at 
MAMI over a wide angular range below pion threshold |^ . It will be interesting to perform 



a best fit for a — (3 and 'y-n- with such an extended data base. 

As one moves to energies above pion threshold, the Compton cross section rises rapidly be- 
cause of the unitarity coupling to the much stronger pion photoproduction channel. There- 
fore this higher energy region is usually considered less 'pure' to extract polarizabilities 
because the procedure would require a rather precise knowledge of pion photoproduction. 
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With the quite accurate pion photoproduction data on the proton that have become avail- 
able in recent years, the energy region above pion threshold for the Compton channel could 
however serve as a valuable complement to determine the polarizabilities, provided one can 
minimize the model uncertainties in the dispersion formalism. In this work, we use the most 
recent information on the pion photoproduction channel by taking the HDT multipoles 
at energies < 500 MeV and the SAID-SP98K solution ||2^ at higher energies. In addition, 
as previously shown in Fig. ^, the subtracted dispersion relations are practically saturated 
by the one-pion channel for photon energies through the A region, which minimizes the un- 
certainty due to the modeling of the two-pion photoproduction channels. In Figs. and O 
we display the sensitivity of the Compton cross sections to 7,^ and a — /? in the lower part of 
the A region in comparison with the available data. As can be seen from Fig. these data 
are quite sensitive to the backward spin polarizability 7,r- This sensitivity was exploited in 
Ref. within the context of an unsubtracted DR formalism, and the value — —27 was 
extracted from the LEGS 97 data, which are shown at the higher energies in Fig. |T^. Our 
results for the subtracted DR are obtained in Fig. [1^ by variation of 77r at fixed a — (3 = 
10, and by variation of a — /3 at fixed 77r = -37 in Fig. |T3|. For 77r we show the results for 
values between 7^ = —27 and 7,^ = —37. One sees that the lower energy data {E^ = 149 
MeV and 182 MeV) can be easily described by the larger values of 7^ if a — /3 decreases to 
some value below 10. On the other hand, the higher energy data [E^ = 230 MeV and 287 
MeV) seem to favor a smaller value of Jt^, and so far we confirm the conclusion reached in 
Ref. P]. However, we have to point out that, at these higher energies, the data around 90° 
cannot be described in our subtracted DR formalism for reasonable values of a — f3 and 7,r- 
This is also seen in Fig. |n| at two fixed angles, now as function of the energy throughout 
the A resonance region, for the MAMI data at 75° and 90° 0. It is again obvious at 
these angles that the sensitivity to 7^ is quite small. Therefore the physics of these data 
is basically driven by pion photoproduction. With the multipoles used here, the 75° data 
are well described, but our prediction falls below the 90° data on the left shoulder of the 
A resonance. In the same energy region there exist also both differential cross section and 
photon asymmetry data from LEGS by use of the laser backscattering technique. In 
Fig. ^ we compare our predictions with these data. One finds that at both energies {E^ 
= 265 MeV and 323 MeV) our subtracted DR formalism provides a good description of 
the asymmetries which display only little sensitivity on 7,^, but underestimates the absolute 
values of the cross sections. In particular close to the resonance position at E^ = 323 MeV, 
the subtracted DR formalism does not allow us to find any reasonable combination of 7^ and 
a — (3 to describe these data. Therefore, within the present subtracted DR formalism, the 
actual data situation at these higher energies does not seem to be conclusive to reliably decide 
on a value of 7,^- Since the uncertainties due to two-pion and heavier meson photoproduction 
are less than 1 % in our subtracted DR formalism, the only possibility to describe the E^ = 
323 MeV LEGS data would be an increase of the HDT Mi+ multipole by about 2.5 % (see 
the dotted hues in Fig. |15D. Indeed such a fit was obtained by Tonnison et al. |§] by use of 
the LEGS pion photoproduction multipole set of Ref. PS| for photon energies between 200 



and 350 MeV and the SAID-SM95 multipole solution ||29[ outside this interval. However the 
more recent SAID-SP98K solution is in very close agreement with the HDT multipoles in 
the A region and hence the predictions also fall below the data at 323 MeV in Fig. p!5| . 
Before coming to any conclusions, we like to point out that new Compton data on the 
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proton in and above the A-resonance region and over a wide angular range have been 
measured recently at MAMI and reported preliminary |^|. These new data will be most 



valuable to check if a systematic and consistent trend becomes visible between the data sets 
at the lower energies and in the A region. 

Finally, in Fig. |16] we show that double polarization observables will be ultimately neces- 
sary in order to reliably extract the polarizabilities a— (3 and 7,^. In particular, an experiment 
with a circularly polarized photon and a polarized proton target displays quite some sen- 
sitivity on the backward spin polarizability 7^, especially at the somewhat higher energies 
of ^ 230 MeV. Such a measurement is more selective to 7,^ due to a lesser sensitivity 
to a — /3 (see Fig. |I6]). It will be indeed a prerequisite to disentangle the scalar and vector 
polarizabilities of the nucleon. 



VI. CONCLUSIONS 

In this work we have presented a formalism of fixed-t subtracted dispersion relations for 
Compton scattering off the nucleon at energies < 500 MeV. Due to the subtraction, the 
s-channel dispersion integrals converge very fast and are practically saturated by the vrA^ in- 
termediate states. Because of the use of subtracted DR we have minimized the uncertainty 
from multi-pion and heavier meson intermediate states. Hence this formalism provides a 
direct cross-check between Compton scattering and one-pion photoproduction. We have 
described this dominant 7A^ ttN —>■ 'jN contribution by using the recent pion photopro- 
duction multipoles of HDT. 

To calculate the functional dependence of the subtraction functions on the momentum 
transfer t, we have included experimental information on the t-channel process through tttt 
intermediate states as 77 — > tttt — >• NN. We have constructed a unitarized amplitude for 
the 77 — > TTTT subprocess and found a good description of available data. In this way, we 
have largely avoided the uncertainties in Compton scattering associated with the two-pion 
continuum in the t-channel, which was often modeled through the exchange of a somewhat 
fictitious cr-meson. 

As the polarizabilities directly enter as subtraction constants in the subtracted DR for- 
malism, it can be used to extract the nucleon polarizabilities from data with a minimum 
of model dependence. We have demonstrated the sensitivity to the polarizabilities a — (3 
and 77r of existing Compton data on the proton both below pion threshold as well as in 
the A resonance region. The effects of the polarizabilities a — (3 and 7,^ on the Compton 
cross sections are strongly correlated. Hence these polarizabilities can only be determined 
simultaneously from unpolarized observables even below pion threshold. When comparing 
the subtracted DR formalism with the existing data in the A resonance region, we have 
found that the actual data situation at these higher energies does not seem to be conclusive 
to reliably decide on a value for 7^. However new Compton data on the proton both below 
pion threshold and in the A region are actually under analysis and will extend considerably 
the experimental data base to fit the proton polarizabilities. 

We have also argued that double polarization observables will ultimately be necessary to 
extract reliably a — (3 and 77r. In particular we have shown that experiments of polarized 
photons on polarized protons show rather large sensitivity to '-/t^ at energies around E^ ^ 
230 MeV and at backward angles. Therefore such polarization experiments hold the promise 
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to disentangle scalar and vector polarizabilities of the nucleon and to quantify the nucleon 
spin response in an external electromagnetic field. 
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APPENDIX A: THE MANDELSTAM PLANE PHYSICAL AND SPECTRAL 
REGIONS FOR COMPTON SCATTERING 



The kinematics of Compton scattering, ■y{q)N{p) l{(l')N{p'), can be described in terms 
of the famihar Mandelstam variables, 

s={q + p)\ t = {q~q')\ u = {q ~ p')\ (Al) 

with the constraint 

s + t + u = 2M^ . (A2) 

Furthermore we introduce the coordinate u perpendicular to t, 

s — u ^ t , , , 

u = — — = E^ + — . A3 

AM ^ AM ^ ^ 

In these equations, is the photon energy in the lab frame and M the nucleon mass. 

The boundaries of the physical regions in the s, u and t channels are determined by the 
zeros of the Kibble function $, 

^s,t,u) =t{us - M^) = . (A4) 

The 3 physical regions are shown by the horizontally hatched areas in Fig. 1. The vertically 
hatched areas are the regions of non-vanishing double spectral functions. These spectral 
regions are those regions in the Mandelstam plane where two of the three variables s, t and 
u take on values that correspond with a physical (i.e. on-shell) intermediate state. The 
boundaries of these regions follow from unitarity. As discussed in Ref. [^, it is sufficient 
to consider two-particle intermediate states in all channels. Since these boundaries depend 
only on the masses, they are the same for all 6 amplitudes Ai. In the Mandelstam diagram 
of Fig. 1 they are symmetric to the line z/ = due to crossing symmetry. For the spectral 
function psu we obtain the boundary 

bi{u,s) = bi{s,u) = [s- {M + m^f][u - {M + m^f] - {ml + 2Mm^f = , (A5) 

and for the spectral function pst we find 

bjj{s, t) = {t- Ami) [s-iM + m,)2] [s - (M - m^f] - 8ml{s + M^ - ml/2) = (A6) 

The boundary of the spectral function put follows from crossing symmetry. We also note that 
these boundaries are obtained for the isovector photon, which couples to a tt+tt^ pair. The 
corresponding boundaries for the isoscalar photon are inside the boundaries of Eqs. ( |A5|) 
and (|A6|) , because it couples to 3 pions. 

APPENDIX B: T-CHANNEL HELICITY AMPLITUDES FOR COMPTON 

SCATTERING 

The t-channel helicity amplitudes for Compton scattering can be expressed in the orthog- 
onal basis of Prange [EO] in the following form: 
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X u{v\\n) \ ^ [Ti + 



P' 

iV2 



+ 



pi^N" - P' m 

+ « 75 ^5 

P' 

p'^N'^ + p'^Nf^ 
+ ^ 75 ^Te ;> vi-p, X^) , (Bl) 



where 



= e^^^^P'^Qpk, , Q = ^{-p-p') = ^{-q' -q), (B2) 

and using the convention 60123 = +1- 

In the cm. system of the t-channel process 77 NN (see Fig. ^ for the kinematics), we 
choose the photon momentum qt (hehcity X'^) to point in the ^-direction and the nucleon 
momentum p' = pt in the xz plane at an angle 6t with respect to the 2;-axis (the anti-nucleon 
momentum is then given by —p = —pt)- In this frame, the helicity amplitudes of Eq. (|BlD 
can be cast into the form 

- ^ (A'^ + A^) 75 T5 



l(A%-A^)75|gl|7'T6| v{-pt,Xf,). (B3) 



2 

Under parity transformation, these amplitudes behave as 

TU^,x',x, M = (-l)A^^-A, T!,^_^^ (^,t) , (B4) 

with the helicity differences and A^r given by = X'^ — X^ (with A^ = or 2) and 
Ajy = X^ — Xj^ (with Atv = or 1) respectively. 

However, the invariant amplitudes Tj(i = 1,...,6) of Prange have the disadvantage to 
behave differently under s ^ u crossing. While Ti, T3, T5 and Tg are even functions of u, 
T2 and T4 are odd functions (note that u — > —u is equivalent to s m). Therefore, L'vov 
used a new set of invariant amplitudes Ai{i = 1, 6), which are all even functions of u and 
at the same time free of kinematical singularities 
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A, = - [Ti + Ts + u{T, + T4)] 



A2 
A3 

A, 
A^ 



t 



[2T5 + Z/(T2 + T4)] 



M2 



M4 - SM 

1 



Ti - T3 



4z/ 



(T2 - T4 



2MT6 - — (Ts 
4i/ 



4. '^^ + ^' 



In terms of the t-channel helicity amplitudes T^^^^ y 
plitudes Ai{u,t) {i = 1,...,6) of Eq. (^) read 



(B5) 

V, t) of Eq. ( pi| ) the invariant am- 
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t Vt - 4M2 
1 



11 11 + -'ii -1-1 

9.''- ^ 9.9'-^ ^ 
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A, 
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(B6) 



In the subtracted DR of Eq. (]T6|), the t-channel integral runs along the line u = 0. There- 
fore, we have to determine the imaginary parts ImtAilu = 0,t) of the invariant amplitudes 
of Eq. (|B(j|). We start by decomposing of the t-channel helicity amplitudes for 77 NN 
into a partial wave series, 



2J+ 1 



(B7) 



where dj^^j^^ are Wigner (i-functions and is the scattering angle in the t-channel, which is 
related to the invariants u and t by 



COS^^f 



4Mz/ 



ViVt- 4M2 ■ 



[B8) 
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It is obvious from this equation that z/ = corresponds to 90° scattering for the t-channel 



process. As explained in Section we calculate the imaginary parts of the t-channel helicity 



amplitudes ^AjvAj^.a' a^(^; ^) through the unitarity equation by inserting vrvr intermediate 
states, which should give the dominant contribution below KK threshold, 

1 l^^l 



[4:71) \Jt J 

Combining the partial wave expansion for 77 — > tttt , 



■KIT NN . 



(B9) 



rp77 — > TTTT 



(t, e^^) = ^ — ^ t; 

Jeven 



J (77 — >■ tttt) 



it) 



(J + A,)! 



(BIO) 



and the partial wave expansion for vrvr — ^ A^A^, 



(«) 



1^ ■ ■ 



We can now construct the imaginary parts of the Compton t-channel partial waves, 



2ImT 



AjvAjviA!yA^ V / 



The partial wave amplitudes Tj 



1 Pit 

Sir) 

An 



T, 



J (77 tttt) 



it) 



T. 



J (tvtt^NN) 



A. 



it) 



(Bll) 



(B12) 



of Eq. (pil| ) are related to the amplitudes f±it) 



of Frazer and Fulco by the relations 



J (tttt^NN) 



Ajv=0 



it) 



IQtt 



ivNP.)' ■ fiit) , 



Pn 



T. 



J (tttt NN) 



Vt 



Aw=l 



it) = 87r^{pNP^)'' ■ fi{t) 
Pn 



with and the cm. momenta of nucleon and pion respectively, 

Pk 



PN = W4 - M2 



9, = ^/t/A 



(B13) 



(B14) 



For the reaction 77 tttt, we will use the partial wave amplitudes Fj\^(t), which are 
related to those of Eq. ( |B10|) by 

2 



T, 



J (77 — > tttt) 



it) 



V2J+1 



■FjA.it). 



(B15) 



Denoting the Born partial wave amplitudes for 77 tt^tt by Bjj>^^{t), the lowest Born 
partial waves (s and d waves) are 



Booit) = 2e' 



1-/3' 



In 



1 + /3 
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1-/32 [ 3-/32 ^1 + ^ 

4 



52o(t) = 2e2 ^ — <! In 



/32 



/5 



1-/3 



4^2 



In 



1-/3 



2 



(BI6) 
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with (3 = p^l {Vi/2) the pion velocity. 

Inserting the partial- wave expansion of Eq. (p^ ) into Eq. (p6|), we can finally express the 
27r t-channel contributions Im(y4j(z/ = 0, t)^'^ by the partial wave amplitudes for the reactions 
77 —>■ 7171 and tttt — *• NN, 



lmtAi{u = 0,t) 



2n 



Pn_ _J_ 



J2 {p.PnY V2jTlFj^_^=,{t) fi*{t) 



J=0,2,4,... 



^ ' J\\ 



ImtA2(z/ = 0,t)2" = 0, 
ImtA3(z/ = 0,t)2- = 

^,T^ (p^PnY V2JTT FjA_,=,{t) 



J=2,4,. 



(_l)('^-2)/2 



(J+1)J (J-i: 



(J-l)(J + 2) J!! 
(J + 2)(J-1) - 2 



(_l)(./-2)/2 2 (j_2)(j + 3) (J - 1)!! 



ImtA5(i/ = 0,t)2- = 
^ /\/f 

ImtA6(i/ = 0,t)2- = 

- — ^ E (P-P^)' ^^27+1 F^A,=2(t) 



V(J + 2)(J-1) 

(J+1)!! 
v/J(J + l) (^-2)!! 



J!! 



J=2,4,, 



X 



(_l)(J-2)/2 [(j + 2)(J-l) - 2] (J-1)! 



v/(J + 2)(J-l) 



J!! 



(B17) 



We note that the s-wave (J = 0) component of the 27r intermediate states contributes only to 
Ai. The amplitude ^2, corresponding to the exchange of pseudoscalar mesons (dominantly 
7r°) in the t-channel, gets no contribution from 27r states, because the 27r system cannot 
couple to the nucleon through a pseudoscalar operator. Furthermore, it is found that only 
waves with J > 4 contribute to the amplitude A^. In our calculations we saturate the 
t-channel dispersion integral with s(J = 0)- and d(J = 2)-waves, for which the expressions 
of Eq. (|B17| ) reduce to those given in Eq. (|^) . 



APPENDIX C: F2 -MESON CONTRIBUTION TO THE 77 ^ tttt PROCESS 

A particle with mass m and spin-2 is described in terms of a symmetric and traceless field 
tensor satisfying the Klein-Gordon equation. Furthermore, as for a massive spin-1 field, 
the 'Lorentz gauge' condition requires that the four- divergence with respect to one of the 
four-vector indices vanishes, 
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( □ + ) ^^'^ = , 
= , 

•l)^'^ = d^^f''' = . (CI) 

Therefore, the tensor ^'^'^ has only five independent components. 

A state of spin-2 is characterized by its polarization tensor e^'^{p,A), where A = 
{2, 1, 0, —1, —2} defines the polarization, 

e^'^ {p, A) = e^^ {p, A) , 
£%(p,A)= 0, 

p,e^''ip,A)= p,e^''ip,A) = 0. (C2) 

The polarization tensor e'^" {p, A) can be constructed from (massive) spin-1 polarization 
vectors by 

£'^^(p,A) = Yl E (lA, lA'|2A)e'^(p,A)e^(p,A'). (C3) 

A=-1,0,1 A'=-1,0,1 

If one chooses the 2;-axis along the momentum of the particle, the three polarization vectors 
of a massive spin-1 particle are 

e^{p,\ = +1) = (0,--^,--J=,0), 

e'^(p,A = -1) = (0,+i=,--^,0), 

e^{p,X = 0) = (b1, 0, 0, /)/m, (C4) 

where p'^ = (p°,0,0, \p\) and p"^ = iv?. In the case of the spin-2 polarization tensor con- 
structed as in Eq. (|(J3D, the following normalization and completeness conditions hold : 

e^'' {p, A) ■ e;, {p, A') = 5aa' , 

J]e'^'^(p,A) [e"^(p,A)]* = - ( K^" K''^ + K'''' K"'^ - -K'^'K"^) , 

with K^" = -g^" + ^ . (C5) 
Finally, the propagator of the spin-2 particle with total width F takes the form 

We will now apply this spin-2 formalism to describe the s-channel exchange of the /2 
meson in the process 77 — tttt. 

The coupling of the (isoscalar) /2(1270) meson to a pion pair (with momenta p^r, p'^ and 
cartesian isospin indices a, h) is described by the amplitude 
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M ih ^ vttt) = ^ 6a, p'/ p: e,. ip, A) , (C7) 

where p is the /2 meson momentum and rrif^ its mass. The couphng constant gf^T^n is 
determined from the /2 —>■ vrvr decay width : 

/2 



where pj^ = ^'w?j^/'i — is the pion momentum in the /2 rest frame. Using the partial 

width r(/2 ^ tttt) = 0.846 To and the total /2-width Tq = 185 MeV l2|], Eq. Q yields 
for the coupling : gj^T,-,, ^ 23.64. 
The Lorentz structure of the vertex /2 77 is given by 

M {h - 77) = -z 2 ^ T^' {q, A,) (g', A;) s,, {p, A) , (C9) 

where JF"^ is the electromagnetic field tensor : 

^"^(g,A,) = g°£^*(g,A,) - g^£°*(g,A,) . (CIO) 
Using the vertex of Eq. (|CDD, the /2 — > 77 decay width is calculated as : 

-77) = ^9%..^h- (cii) 

Using the partial width r(/2 77) = 1.32 ■ 10^^ Fq |2^, Eq. ( pil|) determines the value 
of the coupling constant : (7/277 — 0.239. 

Using these couplings and vertices, we can now calculate the invariant amplitude for the 
process 77 — > /2 tttt : 

Mi^j^nn) = -z2e^^^'^tg,A7mV,A;) A^.„;,(p,A)^p>;^ (C12) 

To determine the 77 mi helicity amplitudes defined in Eq. (|32D, we shall evaluate 
Eq. (|C12|) in the cm. system. For the case of equal photon helicities (A7 = 0) the /2 does 
not contribute, i.e. 

<to = . (C13) 
For the case of opposite photon helicities (A^ = 2) we find after some algebra : 

where ^tttt is the pion cm. angle and f3 the pion velocity as in Eq. (^). It is immediately 
seen from Eq. ( |C14]) that the /2 meson contribution to the d-wave is given by : 

p(/2) (+\ — _J ^ 9f2ii 9f2mT t(3 f^^K\ 

J^j=2A^=2[t) - \/ 15 4 ^2^ t-m%+ tmf.To ' ^ ' 
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TABLES 



TABLE L The contribution of the dispersion integrals to the spin polarizabiUties of the proton. 
The set HDT(l7r) is calculated from the one-pion photoproduction multipoles of the HDT analysis 
1 28 1, while the column HDT gives the total results with the additional contribution of inelastic 
resonance channels. The entries in the last column are the predictions of the dispersion calculation 
of Ref. m. 



7j — excit. 


HDT(l7r) 


HDT 


Ref. g 


(p) 
7i 


+4.83 


+4.33 


+3.1 


(p) 


-0.81 


-0.74 


-0.8 


(p) 

73 


-0.30 


-0.02 


+0.3 


7?^ 


+3.19 


+2.93 


+2.7 




-0.75 


-0.80 


-1.5 


7i? 


+4.23 


+4.29 


+3.7 


7if 


-1.56 


-1.53 


-2.3 


7^^^ 


+10.41 


+9.46 


+7.8 



TABLE II. The same as in Table I in the case of the neutron. 



7j — excit. 



7i 

(n) 
72 

(n) 
73 

(n) 

71 



HDT(l7r) 



HDT 



Ref. |2| 



+7.10 
-0.68 
-1.04 

+3.92 



+7.00 
-0.68 
-0.99 

+3.88 



+6.3 
-0.9 
-0.7 

+3.8 



7o 

(n) 

7l3 
in) 

7i4 
(n) 



-0.06 
+5.02 
-0.74 
+14.27 



-0.09 
+5.02 
-0.77 
+14.09 



-0.4 
+4.9 
-1.3 
+13.0 
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FIGURES 




FIG. 1. The Mandelstam plane for Compton scattering. The physical regions are horizontally 
hatched, whereas the spectral regions are vertically hatched. 
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2 Im 




FIG. 2. t-channel unitarity diagrams for Compton scattering. 
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y + p ^ y + p ( e^iab = 180°, = 250 MeV) 
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FIG. 3. Convergence of the s-channel integral for the amplitudes A\ and ^2- Results for the 
unsubtracted dispersion integral of Eq. (|lO| ) for the one-pion channel (dotted lines) and including 
the two-pion channel (dashed-dotted lines) in comparison with the subtracted dispersion integral 
of Eq. ( |l5|) for the one-pion channel (dashed line) and including the two-pion channel (full lines), 
as function of the upper integration limit z^upper- 
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FIG. 6. Total cross section for the 77 tt^tt" process as function of the cm. energy : 
Born terms (dotted line), Born amplitude with unitarized s-wave (dashed-dotted line), /2(1270) 
resonance contribution (dashed line) and total amplitude (full line). 
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FIG. 7. Differential cross section at various cm. energies for the 77 tt+tt" process : Born 
terms (dotted line), Born amplitude with unitarized s-wave (dashed-dotted line, only shown at the 
four lower energies), and total amplitude including the /2(1270) resonance contribution (full line). 
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FIG. 8. Convergence of the t-channel integral for the amplitude Ai. Results for the unsub- 
tracted (dashed curve) and the subtracted (full curve) t-channel dispersion integrals are shown as 
function of the upper integration limit fupper- Both results are normalized to their respective values 



at tupper = 0.78 GeV^. 
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Y + p^Y+p(t = - 0.163 GeV^) 
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FIG. 9. Real parts (full lines) of the subtracted s-channel integral of Eq. (^) and imaginary 
parts (according to Eq. (|2^)) to the invariant Compton amplitudes Ai,...,Aq as function of u at 
fixed t = -0.163 GeV^. 
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y+p^Y + P (ot-P=10) 




FIG. 10. Differential cross section for Compton scattering off the proton as function of the lab 
photon energy and at 4 scattering angles Q^J^^. The Born result is given by the dotted lines. The 
total results of the subtracted dispersion formalism are shown for fixed a — /3 = 10 and different 
values of "f-K '■ = -37 (dashed-dotted lines), 7^ = -32 (full lines) and 7^ = -27 (dashed lines). 
The data are from Ref. ||l| (circles), Ref. (triangles) and Ref. (squares). 
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FIG. 11. Differential cross section for Compton scattering off the proton as function of the 
lab photon energy and at 4 scattering angles Q^^^ as in Fig. IC. The Born result is given by 
the dotted lines. The total results of the subtracted dispersion formalism are presented for fixed 
77r = —37 and different values ofa — /?: a — /3 = 10 (dashed-dotted lines), a — /3 = 8 (full lines), 
and Q — /? = 6 (dashed lines). Data as described in Fig. Il0|. 
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FIG. 12. Differential cross section for Compton scattering off the proton as function of the 
cm. photon angle for different lab energies. The total results of the subtracted DR formalism are 
presented for fixed a — /3 = 10 and different values of ^-k '■ ^-k = -37 (dashed-dotted lines), = -32 
(full lines) and = -27 (dashed lines). The data are from Ref. Q (solid circles), Refs. (open 
circles) and Ref. Q (squares). 
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FIG. 13. Differential cross section for Compton scattering off the proton as function of the 
cm. photon angle for different lab energies as in Fig. The total results of the subtracted DR 
formalism are presented for fixed ^-j^ = —37 and different values of a—/? : a— (3 = 10 (dashed-dotted 
lines), a — /3 = 8 (full lines), and a — /3 = 6 (dashed lines). Data as described in Fig. 
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Y + p ^ Y + p 




FIG. 14. Differential cross sections for Compton scattering off tlie proton at fixed cm. scatter- 
ing angle through the A resonance region. The total results of the subtracted DR formalism are 
shown for fixed a — = W and different values of 'Jtt ■ = -32 (full lines) and = -27 (dashed 
lines). The MAMI data are from Q (upper panel) and [^j (lower panel). 
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FIG. 15. Photon asymmetries (upper panels) and differential cross sections (lower panels) for 
Compton scattering off the proton in the A resonance region. The total results of the subtracted 
DR formalism are shown for fixed a — P = 10 and different values of 'jn ■ = -37 (dashed-dotted 
lines), = -32 (full lines) and = -27 (dashed lines). We also show the result for a — /3 = 10 
and = -32 when increasing the HDT Afi-|- multipole by 2.5 % (dotted lines). The data are from 
LEGS 97 ||. 
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Y + p ^ Y + p 




FIG. 16. Double polarization differential cross sections for Compton scattering off the proton, 
with circularly polarized photon and target proton polarized along the the photon direction (upper 
panels) or perpendicular to the photon direction and in the plane (lower panels). The thick (thin) 
lines correspond to a proton polarization along the positive (negative) direction respectively. The 
results of the dispersion calculation are for a — fi = 10 and different values for 7,^ : = -32 (full 
lines), = -27 (dashed lines) and 7,^ = -37 (dashed-dotted lines). We also show the result for 
a — /? = 8 and = -37 (dotted lines) . 
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